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Abstract 

The Einstein equations with a positive cosmological constant are coupled 
to the pressureless perfect fluid matter in plane symmetry. Under suitable 
restrictions on the initial data, the resulting Einstein-dust system is proved 
to have a global classical solution in the future time direction. Some late 
time asymptotic properties are obtained as well. 



1 Introduction 

In [T3] plane symmetric solutions of the Einstein- Vlasov system with positive cos- 
mological constant were investigated. It was shown that a spacetime of this type 
which is initially expanding exists globally in the future when expressed in an 
areal time coordinate t and information was obtained about its asymptotics for 
t — > oo. It is future geodesically complete and resembles the de Sitter solution at 
late times. Information is obtained on the decay rates of the components of the 
energy-momentum tensor. 

This paper is concerned with the question, to what extent analogues of these 
results for the Einstein- Vlasov system hold in the case of the Einstein-dust system. 
There is an issue which has to be addressed right at the start. This is that of shell- 
crossing singularities. The intuitive idea behind this concept, as explained in [5] and 
[9 , is the following. A shell of dust particles which are related to each other by the 
symmetry of the spacetime moves in a coherent way. If two of these shells collide 
then the intermediate shells arc trapped between them, so that the matter density 
is forced to blow up. For more information about shell-crossing singularities see for 
instance [7] and references therein. For general plane symmetric solutions of the 
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Einstein-dust system it must be expected that shell-crossing singularities develop, 
even from smooth initial data. For that reason global classical solutions cannot be 
expected to exist. In order to have a global existence theorem in the framework 
of classical solutions it is therefore likely to be necessary to make some kind of 
smallness assumption on the initial data. There is another fact which is essential 
in the following. This is the presence of a positive cosmological constant, known 
as a mechanism leading to solutions of the Einstein equations with exponential 
expansion. Results related to this idea have been obtained in [T^] and [TO]. They 
are used in this paper to guess the decay rates on the geometric and matter 
quantities providing a basis to the bootstrap argument for the proof of the main 
result (Theorem 13.31) . 

The rest of the paper is organised as follows. In section 2 the equations are 
derived and a local existence theorem is obtained for the corresponding Cauchy 
problem. In section 3 the solution is shown to exist globally in the future time 
direction and its late-time asymptotic behaviour is investigated, provided some 
restrictions on the initial datum. 

2 Preliminaries 

2.1 The Einstein equations 

Let (M,g) be a spacetime, where the manifold is assumed to be M — I x T 3 , 
I is a real interval and T 3 = S 1 x S 1 x S 1 is the three-torus. The metric g and 
the matter fields are required to be invariant under the action of the Euclidean 
group E% on the universal cover. It is also required that the spacetime has an 
£?2-invariant Cauchy surface of constant areal time. In such conditions the metric 
is assumed to have the form 

ds 2 = _ e 2v(t,x) dt 2 + e 2X(t,x) dx 2 + t 2 {dy 2 + ^ ^ 

where t > 0, r\ and A are periodic in x, and y and z range in [0, 2ir]. 
The Einstein equations read 

G af3 + Kg al3 = 8irT af3 , (2.2) 

where G Q/3 is the Einstein tensor, T af3 the energy-momentum tensor and A is 

the cosmological constant we assume to be positive. We introduce the notation 
p = e 2.) T 00 ; j = e \+ VT oi and s = e 2A T n_ 

After computations in plane symmetry with the previous coordinates consid- 
erations, we obtain from (|2.2[) the following equations where the subscripts t and 
x refer to partial derivatives with respect to t and x respectively: 

e' 2ri {2t\ t + 1) - At 2 = 8irt 2 p (2.3) 

e~ 2 ' ) (2ir ?t - 1) + At 2 = 8nt 2 S (2.4) 
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e 



Vx = -Ante x+r >j (2.5) 
!A (v™ + Vx(Vx - A*)) - e~ 2 " [\ tt + (A* - »b)(A t + ^ + A = 0. (2.6) 



2.2 The equations for dust 

We consider a pressureless perfect fluid with energy density p := /i(t, x) > and 
4- velocity U a . The latter is normalized to be of unit length U a U a = — 1. The plane 
symmetry allows us to set U a :— £,(e~ v , e _A u, 0, 0) where £ = (1 — u 2 )^ 1 / 2 is the 
relativistic factor, u :— u(t,x) being the scalar velocity that satisfies \u\ < 1. The 
energy momentum tensor for a pressureless perfect fluid is 

T af3 = pU a U> 3 , 

that is 



e-^j 



e- 2X S, 



rpOO 


= e 


-2n 






1 — u 2 ' 




= e7 


X- 


1-u 2 


/Jill 


= e~ 


-2A 


flU 2 




1-u 2 • 



the other components being zero. 

The equations for dust are given by 

V a T aP = 0. (2.7) 

The components \7 a T a2 and V a T a3 vanish identically. Computing the remaining 
two components gives 

(e x p) t + {&j) x = -X t e x S - n x e*j - je x p (2.8) 

(e x j) t + {e r 'S) x = -X t e x j - rj^p - Vj, (2.9) 

and expressing the equations (|2.8p - (12.9|) in terms of the variables p and u gives 

(1 - u 2 )[p t + ue n ~ x p x ] + p[2uu t + (1 + u 2 )e r '- x u x ] 

= -fi(l - u 2 )[(l + u 2 )\ t + 2t~ x + 2u7- lx e r) - x ] (2.10) 
(1 - u 2 )u[p t + ue n ~ X p x ] + p[(l + u 2 )u t + 2ue n ~ X u x ] 

= -n(l - u 2 )[2u(X t + r 1 ) + (1 + u 2 )^" - *]. (2.11) 
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Adding and subtracting (|2.10p and (|2.11[) yields 

(1 - u)\pt + ue"~V] + + e^ X u x ] 

= - M [(l - u 2 )(\ t + Vx eJ'- x ) + 2(1 - u^- 1 } (2.12) 
(1 + u)\p t + ue'^V^l + + e v ~ X u x ] 

= -(1 + - «)A t + 2*" 1 - (1 - u^e""*], (2.13) 

and the linear combinations (|2~T2]l + (f2~Hi)) and (1 + m) ([2~T2]) + (u - 1) ([2"7T5]) lead 
to 

D/x = -fj, [(1 - u 2 )A t - 2r x + e 7 '- x u x ] (2.14) 
Du=(l-u 2 )[-u\ t -e r >- x r} x ], (2.15) 

where the derivative is 

D := d t + ue^- x d x . 

This is called the characteristic derivative and we denote the corresponding integral 
curve by (t, 7). This means that this curve satifies the differential equation 
7t = ue n ~ x , and on the curve D = jjj so that we can for instance rewrite (|2.14|) in 
the form 

^(t, 7 (t)) = -ai [(1 - u 2 )A t - 2T 1 + e"" V] (i, 7 (i)). 

2.3 The Cauchy problem and local existence 

The object of our study is the plane symmetric Eintein-dust system (|2.3[) - (I2.6[) . 
(|2.14p - (|2.15p with unknowns A, r), /1 and u. The initial data are prescribed at some 
time t > 0. To analyze the solutions of these equations, the first step is to obtain 
a local existence theorem. The method, which has been used for instance in [5] for 
the Einstein- Vlasov system in plane symmetry, consists on constructing an iter- 
ation and proving its convergence. In the present investigation we follow another 
approach, which is described in [9] . There are several steps and it is not convenient 
to keep track of the differentiability in the process. For this reason only the case 
of C°° initial data will be treated in this paper. There is a general (without sym- 
metry) local existence theorem for the Einstein-dust system by Choquet-Bruhat 

DP- 

In fact in that reference the author writes the Einstein-dust system (|2.2I) - (|2.7I) in 
the equivalent form 

u v XI v G a p = — 87r/zu Q ,M^V !/ M ly (2-16) 
u"V v up = (2.17) 
V^u")=0 ) (2.18) 
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in order to prove that the equations are hyperbolic in the Leray sense. Then by the 
Leray-Dionne theory [6]- [3] this implies existence and uniqueness of the solution 
to the Cauchy problem for Einstein-dust equations. For more details see also [2J. 
We can apply this general result to the case with symmetry. Consider a plane 
symmetric compact C°° initial data for the Cauchy problem. Then the symmetry 
is inherited by the corresponding solutions. (See [3], section 5.6 for a discussion of 
this.) Next areal coordinates can be introduced in the spacetime. The conclusion 
is a local existence and uniqueness theorem for the plane symmetric Einstein-dust 
system. The solution, like the initial data, is C°°. 

3 Global existence to the future and asymptotics 

This section is concerned with the main result of this paper. We first prove two 
lemmas dealing with bounds on the unknowns as well as all their derivatives. 

In what follows C will denote a positive constant estimating functions that 
are uniformly bounded, and it may change from line to line. 
Let P n (t) := ||5«A(t)|U + + |&(*)lloo 

Lemma 3.1 Consider a C°° plane symmetric solution of Einstein- dust system on 
a time interval [0,T) such that the following estimates hold, where w := e u x : 

I A t - t -1 | < Ct-\ |At 2 e 2j > - 3| < Ct- 3 , \X tx \ < Ct~ A (3.1) 

|»7»| < Cr\ \ Vxx \ < Ct- 3 , M < Ct~ 3 (3.2) 

M < Ct' 1 , |Mx| < Ct~ 3 , \u x \ < Ct' 1 , K| < Cr 2 . (3.3) 

If all derivatives with respect to x of order up to n of the quantities A, \i, w, rj 
and u are bounded then all derivatives with respect to x of order up to n + 1 of the 
same quantities are bounded. 

Proof. Note that from the hypotheses of the lemma it follows that the quantities 
A, pb, w, rj and u, as well as all their first order derivatives are bounded on [0, T). 
By definition w — e- x u x , which implies 

u xx = e x w x + \ x u x , (3.4) 

and from (|2~5|) . 

Vxx = -±Trte v+x j x + r) x (r) x + X x ). (3.5) 

Differentiating u xx and r\ xx n — 1 times with respect to x, the boundedness of 
d x +1 rj and d x +1 u follows immediately from the hypotheses of the lemma. 
The expression for X tx is 

A te = Vxe 2,1 (At + 87rfp) + ATTt Px e 2n . (3.6) 
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Differentiating this n times with respect to x gives a iinear equation for d t (d x +1 X) 
with coefficients which are known to be bounded, except for terms involving d x + p. 
But the latter can be estimated in terms of d x +1 p so that the following holds 



\d t {d n x +1 \)\<C{l + \d n x +1 p\). 
Integrating this in time implies that 

||a? +1 A(t)|L < pr 1 A(t )||c O + C f\l + P n+1 (s)) ds. 

Now recall that 

Dfi x = -fi x [(l - u 2 )\ t - 2t~ l + 2e"- x u x + u(r] x - \ x )e^ x ] 
- p[-2uu x X t + (1 - u 2 )X tx + (rj x - \ x )e n ^ x u x + e' n ~ x u xx 

and 



(3.7) 



(3.8) 



(3.9) 



D{w x ) = -ktu{\ - u 2 )(j] xx + 2 V 2 X - \ x r) x )e 2 ^ x - 2r) x we r >[w + tA(l - 3u 2 )e TI ] 
+ 3t~ 1 uu x w(At 2 e 2v - 1) - uw[rj xx + rj x (r) x - \ x )]e v ~ x - X tx w 



A t + [2u x + u{2r] x ~ \ x )]e*>- x + 



fr l {\ - iw 



^(At 2 e 2r > - 1) 



(3.10) 



Differentiating (|3.9p n times with respect to x and using (|3.5[) and (|3.6p shows that 
D(d x +1 p) depends linearly on <9™ +1 A, d x +1 p, and d" +2 u with bounded coefficients. 
But <9" +2 it can be estimated in terms of d™ +1 w and <9™ +1 A. It then follows that 



|D(^+V)I < C(l + l^ + Vl + l^ l+1 A| + \d^w\), 
and integrating this along the characteristic 7 implies that 

||5; i+ VWIIoo < ||^+V(i )IU + c [\i + p n+1 (s)) d s . 

Jto 

Likewise taking the x-derivative n times in (|3.10p leads to 

\D(d: +1 w)\ < C(l + |^ + V| + l^ +1 A| + \d« +1 w\), 
and integration along 7 implies 

R n+ M*)Hoo < ||5; i+1 W (to)||oo+C [ (l+P n+1 ( S )) ds. 



Putting 



(|3"7l2")) and (|3~H| together implies 



fn+i(t) < iVu(to) + C / (1 + ^«+i(s))rfs. 

Jt 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



By Gronwall's inequality it follows that P n +i is bounded and thus so are d x +1 X, 
d x +1 p and d x +1 w. This completes the proof of the lemma. □ 



G 



Lemma 3.2 If the hypotheses of Lemma \3.1\ are satisfied and if all derivatives 
of the quantities X, r\, fi and u of the form <9 t fe <9™ with n arbitrary and k < m 
are bounded then the derivatives of the form d™ +1 d x of the same quantities are 
bounded 

Proof. From the evolution equations we have 

A t = i(Ate 2 " - r 1 ) + 47rfe 2 " ( o (3.16) 

T) t = ^{r 1 - Me 271 ) + 4 7 rfe 2 "S' (3.17) 

fit = -ue"-V* - M [(1 - u 2 )\ t - ir 1 + e^ x u x ] (3.18) 
u t = -ue^ x u x + (1 - u 2 ) [-uX t - e"-V] • (3.19) 

Differentiating (|3.16[) - (|3.1 9[) n times with respect to x and m times with respect 
to t allows <9™ +1 <9™A, d t m+1 <9™?7, <9™ +1 <9> and <9™ +1 <9> to be bounded. □ 
We can now prove the main result of the present investigation. 

Theorem 3.3 Consider any C°° solution of Einstein-dust system with positive 
cosmological constant in plane symmetry written in areal coordinates with C°° 
initial data. Let 8 be a positive constant and suppose the following inequalities 
hold: 

|A t (*o) -^o 1 ! < \At 2 e 2 ^ -3\<6, \n x {t )\ < 5, \X tx (t )\ < 6, (3.20) 

|»?xx(*o)| < M(*o) < l«(*o)l < S, \» x (t )\ < 8, \u x (t )\ < 8, \w x (t )\ < 5. 

(3.21) 

Then if 5 is sufficiently small, the corresponding solution exists on [to,oo). More- 
over, for this solution the following properties hold at late times: 

\X t -^1= 0(i- 4 ), \At 2 e 2r > - 3| = 0(r 3 ), \ Vx \ = 0{r 3 ), (3.22) 

\vxx\ = o(r 3 ), \x tx \ = o{r A ), n = o(r 3 ), (3.23) 
l/xxl = o(r 3 ), \u\ = o{r l ), \u x \ = oit- 1 ), \ Wx \ = o(r 2 ). (3.24) 

Proof. The proof proceeds by a bootstrap argument. 

By continuity it follows from the hypothesis (|3.20[) - (|3.21[) that 

|A*(t) - t- 1 ] < 26, \At 2 e 2 ^ - 3| < 25, \ Vx (t)\ < 26, \ Vxx (t)\ < 26 

\X tx (t)\ < 26, n(t) < 26, \u(t)\ < 26, |/x x | < 26, \u x \ < 26, \w x \ < 26, 
for t close to to. 

Let C\ and e be constants for < Ci < 1 and < e < 1/2. We can reduce 
6 if necessary so that 26 < Cimin(t ~ 4+e , t^ 1+£ ). Then there exists some time 
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interval on which the solution of the Einstein-dust system exists and the following 
bootstrap assumption is satisfied 

I A ( - t~ x \ < di- 4+£ , \At 2 e 2r > - 3| < di- 3+£ , |A tx | < dt- 4+£ (3.25) 
\Vx\ < dr 3 ^, \ Vxx \ < dt- 3+£ , /i < dt~ 3+e (3.26) 
|«| < dt' 1+e , \ti x \ < dt~ 3+e , \u x \ < dt~ 1+£ , \w x \ < dt~ 2+e . (3.27) 

Consider the maximal interval [to,t*) on which the solution of the full system 
IP3 ]) -([2TB ]) . (I2HD-(123SD exists and 1(5755)) - (|5T27)I hold. Suppose t* is finite. 

Putting inequalities (|3.25[) - (|3.27p into equations coming from the system 
(|2.3p - (|2.6p . (|2.f 4p ~ (|2.f 5p allows new estimates to be derived. We first derive an 
estimate for u. For this purpose an evolution equation for tu can be obtained from 
(|2.f 5p . using the field equations (I2.3p - (|2.4p involving At and rj x . The result is 

D{tu) = u 3 - |(1 - u 2 )(Ae 2 "i 2 - 3). (3.28) 

Using the bootstrap assumption on |u| and \M 2 e 2ri — 3|, integrating the resulting 
inequality along the integral curve 7 and keeping the worst powers it follows from 
([3~2H]) that 

\u{t)\ < [*o|t*(*o)| + Cl + CDt' 1 =: C 2 t-\ (3.29) 
Next we derive an estimate for p. Using (12. 14[) an evolution equation for t 3 p follows: 

D(t 3 p) = t 3 p [t^u 2 - (1 - u 2 )(X t - r 1 ) - e'^V] . (3.30) 

An estimate for e 2 '', e 1)_A and e v+ will be also required. 

e 2 V = A -i t -a (At a e 2,j 

< A- 1 <- 2 [(A< 2 e 2 " -3) + 3] 

< ZA^r 2 + C 1 A- 1 r 5+£ . (3.31) 

On the other hand by assumption |A t — t _1 | < C\t~ A+e , and integrating this in 
time implies that e~ A < e^'^hot- 1 and e A < e Cl+x ^ t "H^ 1 t. Thus 

e n-x < A -i/2 ( V3+C 1 1/2 )t e Cl " A(to) t- 2 , and e" +A < A- 1 / 2 (V3+C 1 1/2 )^ 1 e Cl+A ( t ^. 

(3.32) 

Using this, the bootstrap assumption, integration along 7 and keeping the worst 
powers it follows from (|3.30|) that 



It 



< 



4fx + CxC 2 + Cf + CiA" 1/2 (\/3 + C 1 1/2 )toe c ' 1 - A( * o) l t~ 3 =: C 3 i~ 3 . (3.33) 



Now estimates for |Ai 2 e 2,? — 3| and |tAt — 1| will be derived. Estimates for the 
matter quantities S and p are needed for this purpose. From the definition of S 
and p and the estimates for p and u obtained above, we obtain 

P<j^t- 3 , S<^lr\ (3.34) 
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From (|2.4p we have 



d t [--te~ 2v {At 2 e 2r > - 3)] = -87rf 2 S, 

so that using (|3.31[) . (|3.34[) . integration and keeping the worst powers gives 
|Ai 2 e 2,) - 3| 



< [3A- 1 + C1A- 1 ] 



t \Atl - 3e 



-M*o)| 



24nC 3 C 2 
I-C2 



(3.35) 



From 



we have 



t\ t -l = i(Ae 2 "i 2 - 3) + 47rf 2 e 2 "p, 

and using (|3.3ip . (|3.34p . (|3.35[1 and keeping the worst powers yields 

C 3 



\t\ t -M < 



C 4 + 47rA- i (3 + Ci)- 



1-C 2 

An estimate for r\ x will be derived. Recalling that 



i~ 3 =: c 5 r 3 . 



(3.36) 



= -4vrfe A+ "j and .7 = 



fXU 



it follows from (|3~29|) and (j3~32l) - (j3~33l) that 



|t? x | < 4 7 rA- 1 / 2 ((7 1 1 / 2 + V3)e Cl+A ^)-^i 

1 — O2 



= : C 6 t- 3 . 



(3.37) 



The following are estimates for u x and /x x which will be needed in order to get 
better estimates for \ tx and r] xx . 

An evolution equation for tu x can be derived by taking the x-derivative in (|3.28|) : 

D(tu x ) = (3m 2 - l)u x (t\ t - 1) + 3u x ii 2 + t{r) x + \ x )uu x e v ~ x - tu 2 x e v ~ x 

- i(l - u 2 )[u\ tx + {r, x - A^ffee"-* + r ?ra e"- A ] ! (3.38) 

and using the bootstrap assumption, some estimates obtained above as well as 
integration along 7 leads to 



r 



Kl < [*oK(*o)| +4C 2 +4< C 1 A- 1 / 2 (V3 + C 1 1/2 )e Cl - A (*°H2C 1 + \X x (t )\) 

= : C 7 t-\ (3.39) 
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Now taking the ^-derivative in (|2.14[) allows us to obtain an evolution equation 
for t 3 \i x : 

D(t 3 ii x ) = -t 2 {tX t ~ 1)[(1 - u 2 )fi x - 2fiuu x ] - 2t 3 f i x u x e r >- x 

- t 3 (r] x - X x )ufi x e v ~ + t 2 (u 2 fi x + 2fj,uu x ) - i 3 (l - u 2 )fj,X tx 

- t 3 (r, x - \ x )fiu x e^- x - t 3 nu xx eV- x . (3.40) 

So far we know how to estimate all the terms on the right hand side in (|3.40|) 
except the last one in u xx which needs to be worked out carefully. We can obtain 
an evolution equation for u xx by taking the ^-derivative of Du x . This does not 
give a satisfactory result because terms containing X xx occur, and we do not know 
how to estimate them. This difficulty can be overcome if we rather take the x- 
derivative of D(e~ x u x ), which will lead to an estimate for w x (where w = e u x ). 
An estimate for u xx will then be deduced from the following relation which is 
obtained by differentiating the equality w = e~ x u x in x: 

u xx = e x w x + X x u x . (3-41) 

Note that the factor e~ x allows us to eliminate bad terms such as X xx . This device 
has been used in another context [13] . 

The term w x will now be estimated. An evolution equation for w is 
Dw — —Xtw — w 2 e v — ur] x we v ~ x — t^ 1 w(l — 3it 2 ) 

- ^ lw(1 2 ~ 3 " 2) (At 2 e 2 " - 3) - M Vx u{l - u 2 )e 2 ^ x , 
so that differentiating this in x implies 

D(t 2 w x ) = -M 3 u{l - u 2 )^ + 2tjI - X x r] x )e 2v - X - 2t 2 r] x we v [w + tA(l - 3u 2 )e IJ ] 
+ 6tuu x w + 3tuu x w(At 2 e 2r) — 3) — t 2 uw(ri xx + rj x (r] x — X x ))e n ~ x — t 2 X tx w 

- t(tX t -l)-3tu 2 + t 2 [2u x + u(2r lx --X x )}e r >- x + ^ — ^(Ai 2 e 2 "-3) w x . 

(3.42) 

We use the bootstrap assumption, some estimates obtained above and integration 
along 7 to get 

Kl < [t 2 \w x (t )\ + 4C 7 C 1 i e Cl ~ A( ' o) (l + |Ax(*o)|) 
+ 4C 1 2 (1 + A- 1 /2 )(1 + lA^io)!)]^ 2 

=: C s t- 2 . (3.43) 
It then follows from (|3~4Tj) that 

\u xx \ < \c s t 1 e ^+^(*o) + c 7 {d + |A»(to)|)l t- 1 =: O^r 1 . (3.44) 
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This together with (|3.40j) imply that 

Ifel < [*olM*o)l + ic\ 

+ (3d +3C„^ + |A K (i )|)(V3 + C 1 1/2 )A- 1 / 2 toC ie Cl - A ( t «)]r 3 

=: C* 9 ^ 3 . (3.45) 

Estimates for p x and j x are required in order to derive estimates for X tx and rj xx . 
We have 

_ 2uu x p 

px ~ ~ + (i_„2) 2 ; 

using estimates established above gives 

Cg + 2C2C3US 3 
< ^tt; i 



and recalling that 
it follows that 



I A* J < 



(1 - d) 2 
A te = Vxe 2,1 (At + 8-Ktp) + 4irtp x e 271 

Cg + 2C2C3C8 



A- 1 C 6 (3 + C 1 )(A+-^+4 7 r^ 1 



1-Ci 



Ciot 



(3.46) 



Now 



+ pu x 2u 2 U x jJL 



l-U 2 (l-u 2 ) 2 
so that using estimates obtained above gives 



C 2 C 9 + C 3 C 8 (2C 2 + n 



=: d r 



Taking the spatial derivative of (|2.5p gives 

rjxx = -Airte n+ j x + r) x (r) x + X x ). 

It then follows that 

|tfc„| < [AttA- 1 / 2 ^ 2 + V3)e c 'i+ x ^C 3x + C 4 (C 4 + C x + \X x (t a )\) 
=: Cut- 3 . 



(3.47) 



(3.48) 



f-3 



(3.49) 



The constants Ci — G\\ appearing along the proof are all less than or equal to 
C x (g(S) + C 2 ), with C a positive constant and g(6) a positive function of 6 
tending to as S tends to 0. Therefore it is always possible to choose C\ and 
S small enough in such a way that CC\ < 1/2 and Cg(5) < C%/2, and so the 
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constants C2 — C\\ are all less than C\. This together with lemmas [3.1113.21 show 
that all derivatives of A, r\, \i and u are bounded on [0,i*). This means that the 
solution can be extended to a time interval [to,ii) on which (|3.25|) - (|3.27p hold, 
with t\ > i*. This contradicts the maximality of the interval [to,t*). Therefore 
t* = 00 and the proof of the theorem is complete. □ 

Acknowledgement : I am grateful to A. D. Rendall for fruitful suggestions. 
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